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Abstract 

| This is an English translation of the manuscript [1] which was appeared in Surikaiseki 

kjfy Kenkyusho Kokyuroku No. 1055 (1998). The asymptotic efficiency of statistical estimate 

of unknown quantum states is discussed, both in adaptive an collective settings. Aaptive 
| bounds are written in sigle letterized form, and collective bounds are written in limitting 

^ ■ expression. Our arguments clarify mathematical regularity conditions. 

<n : 

1 Introduction 

o : 

■ This is an English translation of the manuscript [1] which appeared in Surikaiseki Kenkyusho 
Kokyuroku No. 1055 (1998). 

| In the estimation of the unknown density operator by use of the experimental data, the 

error can be reduced by the improvement of the desingn of the experiment. Therefore, it is 
natural to ask what is the limit of the improvement. To answer the question, Helstrom[3] 
founded the quantum estimation theory, in analogy with classical estimation theory (in 
the manuscript, we reffer to statistical estimation theory of probability distribution as 
'classical estimation theory'). Often, for simplisity, it is assumed that a state belongs 
to a family Ai = {pe\9 £ © C R m } of states, which is called model and that the finite 
dimensional parameter 9 is to be estimated statistically. 

He considered the quantum analogue of Cramer- Rao inequality, which gives the lower 
bound of mean square error of locally unbiased estimate. This bound, however, is not 
^ ' achievable at all, when the number of the data is finite. 

■ Let us assume that the number n of the data tends to infinite. Then, if some regularity 
conditions are assumed, it is concluded that if the estimate is consistent, i.e., the estimate 
converges to the true value of parameter, the first order asymptotic term of mean square 
error satisfies the Cramer-Rao inequality, and that the bound is achieved for all 9 € 0. 
This kind of discussion is called first order asymptotic theory. 

The quatntum version of first order asymptotic theory is started by H. Nagaoka[5] [6]. 
He defined, in our terminology, the quasi-quantum Cramer-Rao type bound, and pointed 
out that the bound is achieved asymptotically and globally. The proof of achivability, 
however, is only roughly sketched in his paper. In this manuscript, the proof of the 
achievability of the bound is fully written out, and the regularity conditions fo the achiev- 
ability is revealed. In addition, we defined another bound, the quantum Cramer-Rao 
type bound, and showed that the new bound is also achievable, if the use of quantum 
correlation between samples are allowed. 
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2 Preliminaries 



An estimate 9 is obtained as a function 9(u) of data lo G ft to R m . The purpose of the 
theory is to obtain the best estimate and its accuracy. The optimization is done by the 
appropriate choice of the measuring apparatus and the function 9{u) from data to the 
estimate. 

Let <j(R m ) be a a- field in the space R m . Whatever apparatus is used, the data 00 G ft 
lie in a measurable subset B G a(H m ) of Q writes 

Pr{w € B\9} = tr p(9)M(B), (1) 

when the true value of the parameter is 9. Here, M, which is called positive operator 
valued measure (, POM, in short), is a mapping from subsets B C ft to non-negative 
Hermitian operators in H, such that 

MO) = 0,M(tt) = I, 

00 00 

M((J^),= E M (^) (B i nB j = cf>,i^j), (2) 
i=i i=i 

(see Ref.[3],p.53 and Ref.[4],p.50.). Conversely, some apparatus corresponds to any POM 
M [10] [9]. Therefore, we refer to the measurement which is controled by the POM M as 
'measurement M\ A pair (9,M) is called an estimator. 

The classical Fisher information matrix Jq 1 by the POM M is defined, as in the 
classical estimation theory, 



•Jq 



where di = d/d9\ P^f(B) := tr pgM(B), and v is some underlying measure (in the 
manuscript, we assume that for any POM M, there is a measure v in fl such that Pg f -< v 
for all G O). Denote the mean square error matrix of (0M) by Vg[9,M], and, as 
the measure of accuracy, let us take TrGVg[M], where G is nonnegative symmetric real 
matrix. If G = diag(gi, • • • .g m ), TrGYg[M] is weighed sum of mean square error of the 
estimate 9 l of each component 9 % of the parameter. 

Let us define locally unbiased estimator (9, M) at 9 by, 

Eg[6 n , M] :=y>Htr^M(cM = P, (j = 1,- ■ ■ ,m). (3) 
9\uj)tTd kPe M{Au) = S{, (j,k = l,---,m). (4) 



/ 



Then, is caracterized by, 

jM' 1 = mi{Vg[9, M]\9: (9, M) is locally unbiased}, 

and the quasi-quantum Cramer- Rao type bound 
Cg(G) is defined by, 

Cg{G) := inf{TrGV 9 [^, M] \ M is locally unbiased} 
= in^TrGJ^- 1 | M is a POM in H}. 
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Nagaoka pointed out that the quasi-quantum Cramer- Rao type bound is achievable asym- 
totically for every 9 £ 0[5]. Cg(G) is calculated explicitely for several special cases [7] [2]. 

Suppose n-i.i.d. pairs pf n of the unknown state pg are given. The the sequence 
{(§, M n )}, where M n is a POM in H® n , is said to be MSE consistent if the estimate 9 n 
by converges to the true value of the parameter in the mean, i.e., lin^^oo Vg[(6 n , M n )] = 0. 



3 The quasi-classical Cramer-Rao type bound 
3.1 The lower bound 

Let M^), M( n ) be a sequence of the POMs in H., and apply the measurement Mr\\ 
to the first sample, and the measurement M( 2 ) to the second sample, and so on. The 
choice of M(k) is dependent on the outcome ujk-i = w (fc-i)) of Mn M/ fe _ 1 ). 

To reveal the dependency of Mha on tOk-i, we write Mrk)[cOk-i]- 
Let us define the POM M n in H n which takes value in VL n by, 

„ n 

M n (B)= I (g)M (fc) [w fc _i](<L; (fc) ). 

Then the data lj„ is controled by the probability distribution P B Mn {B) = tr pf n M n {B). 
The estimator is said to be asymptotically unbiased if 

(B n )*= (B 9 (0 n ,M)y := ^(^H-^Pf(d^)^O asn^oo, (5) 

(A,)j- = (^(e,M)) := w El[9 n , M n ] ^ 6} as n oo. (6) 

The MSE consistent estimator satisfies (5) always. Therefore, if apropriate regularity 
condisions are assumed so that the defferential, the integral and the trace commute with 
each other, then (6) is also satisfied, and the estimator will be asymptotically unbiased. 

Theorem 1 If {(6 n , M n )} is MSE consistent, and lim n _^ oc , nVg[(0 n , M n )\ exists, 

lim nTrGV 9 [{e n , M n )\ > C e (G), (7) 

n— »oo 

Proof In the almost same manner as classical estimation theory, (6) leads to, 

nV [O n , M n ] > nA n (jf")"' 1 A n (8) 
Elementary calculation leads to, 




where MJ} G M is a POM in TL wich id defined by 




(8) and (9) yield 

Tr GnY e [e n , M n ] > Tr GA n [J e •) 1 A n > C e ( A n GA n ). (10) 
Passing both sides of (10) to the limit n — ► oo, we have the theorem. □ 
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3.2 Estimator which achieves the bound 

The estimator defined in the following achives the equality in the inequality (7) if the reg- 
ularity conditions (B.l-4) are satisfied. The proof will be presented later in the subsection 
3.4. 

First, apply the measurment Mq to y/n samples of unknown state pg, and 
calculate 9 n which satisfies (12). Second, apply the measurement to the 
remaining n — ^fn samples, where Mq is defined by 

Tr G (J^y 1 <C e (G)+e f , (11) 

(11) is satisfied. Then, 9 n is defined to be 9 n (9 n ), where 9 n {9') is defined by, 

n dP M *' 
9 n {9') = argmax V log e (u k ). 

3.3 Regularity conditions 

(B.l) There is a POM M and 9 n which satisfies 

lim Pf"{||0 - 9 n \\ > 5} = 0, V<5 > 0. (12) 

n— >oo 

(B.2) K := supg 6 Q ||#|| is finite. 

(B.3) 9 n {9') achieves the equality in classical asymptotic Cramer- Rao inequality of the 
familily {P e \ 9 £ 0} of probability distributions. 

(B.4) The higher order term of mean square error of 9 n (9') is uniformly bounded when 
||#' — 6*|| < 5\ for some 5 > 0. In other words, for any e > 0,9 G 6, there exists a 
positive real number 5\ > and a nutural number N such that, 



(n - y/n) Tr GY e ,n -TiG (j™ 6 ) 



<e + e, Vn>iV, yd s.t. \\9-9\\ < (5i,(13) 



where \e,n is the conditional mean square error matrix of 9 n when 9 n is given. 
(B.5) For any e > 0, 9 G 6, there exists 82 > 0, such that, 



Tr G (J^y 1 -C e (G) 



<e + e', V0, si. 11^-^11 < £ 2 . (14) 



(B.l) is satisfied almost always, and (B.2) is not restrictive . For 9 n {9') is the maximum 
likelihood estimator of the family {P e of probability distributions, (B.3) is satisfied in 
usual cases. The validity of (B.4), however, is hard to verify. Therefore, in the future, 
this condition needs to be replaced by other conditions. Obviously, (B.5) reduces to the 
following (B.5. 1-2), both of which are natural. 

(B.5.1) The map 9 ^ Cq{G) is continuous. 

(B.5. 2) For any 9', the map 9 [J^ 9 '}^ 1 is continuous. 
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3.4 Proof of achivability 

Theorem 2 If the model M. satisfy conditions (B.l-5) in the following, then we have, 

lim nTiGV e [9 n , M n ] = Cg(G), Vfl G 6, (15) 

n— >oo 

Proof Let us choose <5i,<52 and N so that (13 — 14) are satisfied, and define 5' := 
min(5i, ^2). Then, if n > N, we have, 

nTrGV e [§ n , M n ] 
= nj TrGV e [6 n (e n ), M n ]Pf»(dw) 

TiGV e [e n (e n ), M n ]P^(dLo)+K 2 TvG [ Pf n (du;) 

J\\e-d„\\>5' 



n 


/ 
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< - _ y= I _ [TrGU e 9n \ +e + e' \Pf n {dLo) + nK 2 TiGVf n {\\e -9 n \\ > 5} 

[ (Ce(G) + 2e + e')Pf l ( dto) + nK 2 Tr GP^{\\6 - 9 n \\ > 5} 

n - yjn J\\e-e n \\<s 

< 1 ^{Ce{G) + 2e + e')+nK 2 TTGPf n {\\6 -9 n \\ > 5}. 

(B.l) implies that the third term of last end of the equation tends to n as n — ► 00. 
Therefore, we have, for every e' > and for every e > 0, 

lim nTrGV e [e n , M n ] < Cg(G) + 2e + e'. 

n— >oo 

which leads to the theorem. □ 



4 Use of quantum correlation 

In this section, we consider the minimization of asymptotic mean square error where M n 
runs every POM which satisfies MSE consistensy. Physically, this means we allow the use 
of interactions between sapmles. 

So far, we considerd POM which takes value in S7, or the totality of all the possible 
data. Instead, in this section, we consider POM with the values in R d , for if M is a POM 
with values in S7, M o 6 1 is POM with the values in R d . MSE consistensy is defined in 
the same way as the precedent sections. 

Let Cg(G) denote the quasi-quantum Cramer- Rao type bound of the family {pf n \9 G 
0} of density operators in H® n . Then, the quantum Cramer- Rao type bound C®(G) is 
defined by, 

C e °(G) := lim nC r e \G). 

For Cg(G) > nC%{G) holds true, we have, 

C e (G) > <#(G). 



5 



Theorem 3 // the sequence {M n }^ =1 is MSE consistent, we have, 

lim ntrGYg (M n ) > C®(G). (16) 

Proof In the almost same manner as the proof of theorem 1, we have, 

Ve[M n ] > A n {jr)~ 1 l A n , 

n Tr GV e [M n ] > n tr GA n ( jf » ) " 1 l A n 

> nCZ^AnGAn), 

which approaches (16) as n — > oo. □ 

If the family {pf n \9 G 0} of density operators satisfies (B.l-5), we have the following 
theorem. 

Theorem 4 There is a MSE consistent sequence {M n } of POM such that lim n __ >00 ntr G\Q[M n \ < 
C®(G) + e is satisfied for every e > and for every 6 G 0. 

Proof 

Let us devide n samples into ri2 groups each of which is consist of n\ samples, and 
let MJM, . . . , M?K be a sequence of POMs in H® ni . Apply the measurement MJM to the 

first group pf 1 ™ 1 of samples, and apply to the second samples, and so on. The choice 
of Mn\ is dependent on the outcome of the measurements M^, . . . , M^, 1 _ 1 s. With ni 
fixed, let us approach n2 to oo. Then, theorem 2implies the existence of a MSE consistent 
sequence {M n } of POM which satisfies 

lim ntr GY e [M n ] = lim mn 2 tr GV fl [MJ = mCZHG). (17) 

n— >oo ri2— +oo 17 

For any epsilon, if ni is sufficiently large, limn^oo n tr G\q (M n ) < C®{G) + e is satisfied, 
and we have the theorem. □ 
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